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iy~j ' Abstract 

This paper presents empirical evidence supporting Goldfeld's conjecture on the average analytic 
rank of a family of quadratic twists of a fixed elliptic curve in the function field setting. In 

■ particular, we consider representatives of the four classes of non-isogenous elliptic curves over 
Fg(t) with (q, 6) = 1 possessing two places of multiplicative reduction and one place of additive 
reduction. The case of q = 5 provides the largest data set as well as the most convincing evidence 
that the average analytic rank converges to 1/2, which we also show is a lower bound following 
an argument of Kowalski. The data was generated via explicit computation of the L-function 
of these elliptic curves, and we present the key results necessary to implement an algorithm to 

^ ' efficiently compute the L-function of non-isotrivial elliptic curves over F g (t) by realizing such a 

curve as a quadratic twist of a pullback of a 'versal' elliptic curve. We also provide a reference for 
, our open-source library ELLFF, which provides all the necessary functionality to compute such 

L-functions, and additional data on analytic rank distributions as they pertain to the density 

■ conjecture. 

o 

i ; 

1 Introduction 

•l-H 

Goldfeld's conjecture, in its original form, makes an assertion about a family of elliptic curves over 
a number field and some form of rank. For example, if we fix an elliptic curve E/Q and consider 
the set of its quadratic twists ordered by (increasing) discriminant of the twisting fields, then the 
conjecture asserts that the average rank of the first n curves tends to the limit 1/2 as n tends to 
infinity. In this paper we will fix the rational function field K = ¥„(t) as our base field and consider 
families of elliptic curves over K for which we can calculate each family member's analytic rank. 
Little theoretical progress has been made when we consider the average rank of an 'increasing' 
sequence of curves, and the genesis of this paper lies in a computational project to generate a rich 
data set of L-functions for studying the (empirical) averages. 
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In the bulk of this paper we focus on some algorithms of increasing complexity for explicitly 
calculating the L-function L(E/K,s) for an elliptic curve E/K; the increasing complexity allows 
for increased speed but at the cost of increasing the disk and memory space requirements. In 
contrast to L-functions over number fields, these //-functions have the remarkable property that 
if we define a new variable T by T = q~ s , then we can represent the L-function as a polynomial 
L(E/K, T) G 1 + T-Z[T]. The analytic rank of E/K is simply the order of vanishing at s = 1 (T = 
1/q), hence the analytic Goldfeld's conjecture requires only a miniscule amount of the information 
present in each L-functionQ However, there are many other questions one can ask about how 
L(E/K,T) varies with E/K, so we hope that our database and algorithms will prove useful for 
others. 

For the calculation of a single L-function, one can use the theory outlined in the first two sections 
of section [2j We follow the usual approach for calculating the L-function by expressing it as an 
Euler product; the terms of the product are indexed by valuations in K, i.e. by monic irreducibles 
in the polynomial ring F q [t] and a point at infinity. The Euler product is infinite, but there is a 
finite collection of Euler factors that completely determine L(E/K, T). To compute the L-function 
as efficiently as possible (in this approach) we want to minimize the size of this collection (e.g. by 
using a functional equation) and the cost of computing a single Euler factor. 

Those who have experience computing quadratic twists know that the cost of computing an Euler 
factor for a twist is much cheaper if one knows the Euler factor for the original curve, the work 
being reduced to computing a Legendre symbol. A similar efficiency emerges when we consider 
pullbacks, i.e. when we replace K with a finite extension L/K: almost all of the Euler factors for 
L(E/L, T) can be cheaply calculated from the Euler factors of L(E/K, T). In fact, for each q, there 
is an elliptic curve Eq over F = F g (j) such that every E/K may be written as the combination of 
a pullback of Eq to K (via the embedding F — > K induced by j \-t j(E)) and a twist (by quadratic 
L/K). In particular, the Euler factors for L(E/K,T) may be cheaply computed from the Euler 
factors of L(Eq/F,T), and the second half of section [2] explains this in detail. 

The upshot is that whenever we construct a new curve E/K as a pullback or twist of another curve 
Eq/F, then the additional cost of computing L(E/ K,T) using a precomputed table of sufficiently 
many Euler factors for L(Eq/F,T) is much cheaper than if we computed L(E/K,T) from scratch. 
We have written a library for calculating L-functions whose core routines use the methods we outline 



in section [2] and recently provided the wrapper routines for it to be used within Sage [S + 10j . The 
library is called ELLFF (for elliptic L-functions over function fields), and we describe its basic 
structure in section [3j We used our code to gather data to empirically study Goldfeld's conjecture 
in the function field setting, and we report a summary of this data and observations in section HJ 



1 The Birch-Swinnerton-Dyer conjecture asserts that the algebraic (Mordell-Weil) and analytic ranks are the same, 
so someone who likes looking for points may want to comb our database for curves of rank, preferably at least two, 
and try their hand at producing points. 
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2 Theoretical Framework 



2.1 Elliptic Curves 

Fix a prime power q relatively prime to 6. Let K = ¥ q (t) be the function field of the curve P 1 /F g , 
Ok = F q [t] be the affine coordinate ring of P 1 — {oo}, Ooo = ¥ q [u] be the affine coordinate ring of 
P 1 — {0}, and u = 1/t G K. We identify the set of closed points | IP" 1 1 = {tt} with the set formed 
by the closed point tt = oo together with the monic irreducibles tt G ¥ q [t] of positive degree, and 
we write ¥ n for the residue field. For tt = oo we identify ¥ w with the quotient field ¥ q [u]/u, and 
otherwise we identify ¥ n with the quotient field ¥ q [t]/iT. 

Let E/K be an elliptic curve. Up to a change of coordinates, we may represent our elliptic curve as 
the projective plane curve given by the affine curve y 2 = x 3 + ax + b, where a,b G K, together with 
the point at infinity. The discriminant A = 4a 3 + 27b 2 and j-invariant j = 6912a 3 / A are rational 
functions in t, i.e. elements of K, and A / 0. 

Up to a change of coordinates (x,y) i-> (x/g 2 ,y/g 3 ) with g G K x , we may assume o, b G Ok and 
deg t (A) is minimal because Ok is a principal ideal domain, and we write A„- = A for tt / oo. There 
is a unique integer e such that the substitution (t,x,y) >->■ (l/u,x/u 2e ,y/u 3e ) yields a model of E 
over ¥ q [u] satisfying similar conditions and we write Aqo for the discriminant of this model. We glue 
the two models together over the annulus P 1 — {0, oo} to form the so-called minimal Weierstrass 
model £ — > P 1 ; it is the identity component of the so-called Neron model of E. 

For each tt, we write £/F„- for the fiber of £ — > P 1 over tt. If tt / oo, then it is the projective 
plane curve given by 'reducing modulo 7r' the model for E over Ok, while for tt = oo we use the 
model for E over fi/F^ is a smooth curve if and only if the image of A n in F„- is non-zero, 

otherwise it has a unique singular point. We write M for the finite subset of tt such that fi/Fjr is 
singular with a node, A for the finite subset such that £/F,r is singular with a cusp, and U for the 
open complement P 1 — M — A; they are the loci of multiplicative, additive, and good reduction 
respectively of £ — > P 1 . We recall that if tt / oo G M U A, then tt G M if and only if the image of 
b in ¥ w is non-zero, and a similar criterion holds if tt = oo G M U A. 

We decompose M into the subset M + of tt for which the slopes of the two branches through the 
node of £/F 7r are rational over F„- and M for the complement M — M + ; they are the loci of split 
and non-split reduction respectively. The following lemma gives a criterion for deciding whether 
any tt / oo G M lies in M + or M~ , while for tt = oo one can use the model for E over ¥ q [u] to 
deduce a similar criterion. 

Lemma 2.1. If tt ^ oo £ M , then tt G M + if and only if the image of 6b in ¥ w is a square. 

Proof. Over F,r our affine model specializes to y 2 = (x — c) 2 (x + 2c) for some c / G F T and 
the node lies at (x,y) = (c, 0). The substitution y = s ■ (x — c) and cancellation of (x — c) 2 leads 
to the curve s 2 = x + 2c, the blow up of our original curve at the node. The slopes of the two 
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branches are the s-coordinates of the points (x, s) = (c, s) which lie on this curve, hence s = ±\/3c. 
In particular, the slopes are rational over Fjr if and only if 3c is a square in F^. In terms of the 
singular model, we see that the image of b in ¥ n is 2c 3 , hence 3c is a square in F^ if and only if 
66 = 3c • (2c) 2 is. □ 

In the first three lines of the following table, produced using table 15.1 of [S], we give criteria for 
determining the type of additive reduction E has over ir G A. In the last line of the table we define 
a constant e n which will be used in the next section. 



Kodaira symbol 


I* 

n 


T* 


II 


IV 


IV* 


II* 


III 


III* 


ord 7r (A 7r ) 


6 + n 


6 


2 


4 


8 


10 


3 


9 


j (mod 7r) 


oo 


oo 














1728 


1728 




-1 


-1 


-1 


-3 


-3 


-1 


-2 


-2 



Table 1: Additive Reduction Information for E/K 



2.2 L-functions 



We keep the notation of the previous section and add the assumption that j is non-constant, i.e. it 
lies in the complement K — ¥ q , and we remark that most of what follows extends to the case where 
A (but not necessarily j) is non-constant. 

For each tt £ \U\ and m > 1, we write F^m for the unique extension of F^ of degree m, £(F 7r m) for 
the set of F^m -rational points of E./¥ n , and 



m deg(-7r) 



+ 1 - #£(F W ™). 



For each tt G \U\, the zeta function of £/¥ n is given by the exponential generating series 

(OO rpm\ 
#£(F. m )— . 
m=l / 

It is a rational function in Q(T) with denominator (1 — T)(l — g de §( 7r )T) and numerator 

L(T, 8,/¥ T ) = 1 - a n iT + q dc ^T 2 . 



(2.1) 



Because we assumed j is non-constant, the L-function L(T,E/K) is a polynomial in Z[T] with 
constant coefficient 1 (cf. bottom of page 11 of |KJ) and satisfies 



deg(L(T, E/K)) = deg(M) + 2 deg(^) - 4. 
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It has an Euler product expansion 

l(t,e/k)= y[ l (r dc ^ 7T \s./¥ n y l ■ yi (l-r^^y 1 ■ yi (i+T dc ^y\ (2.2) 

Tre\u\ 7reM+ -n-eM- 

Using (|2.ip and the formal identity 1/(1 — aT) = exp(^^L 1 (aT) n /n) it is easy to show that 



L(r dc gW)£/F ^ )=exp ^ deg(Tr)-^ 



n/ dcg(7r) 

n 

n>l,deg(7r)|n 



Therefore, if we define 



b n = Y, deg(7r)-a WdcgW + Yl deg(^)+ £ deg(7r)(-ir/ de sM 

7rg|C/|,deg(7r)|ri 7reAf+.deg(7r)|n tt<=M-, deg(7r)|n 



then we can rewrite (|2.2|) as 

(oo ™ \ 
)• (2-3) 

If we truncate the formal series expansion of the right side of (|2.3p by reducing modulo T N+1 for 
N > 0, then we obtain the congruence 

L(T,E/K) = exp (E> — J ( mod ^ 

Thus L(T, S/JT) (mod T N+l ) is completely determined by {b n : 1 < n < N} for iV = deg(L(T, E/K )), 
and by definition this set is determined by the Euler factors over tt E |P x | satisfying deg(7r) < N. 
In fact, by taking the functional equation into consideration, as described below, it suffices to take 
N=[deg(L(T,E/K))/2\. 

If we write L{T,E/K) = J2n=o Cn ^ n ^ or ^ = deg(L(T, then to recover cq, ...,cjy from 



(|2.4p it suffices to apply the following lemma. 

Lemma 2.2. If {cq = 1} U {b n ,c n : 1 < n < iV} are numbers satisfying 



N ™\ A 



cxp 



E b «— hE c « T " ( modT 



n , 

vn=l / n=0 



i/ien t/iey satisfy the recursive relation 

if 

n — ^ 



n 

m=l 



Proof. If we take the (formal) logarithmic derivative of both sides of the assumed relation between 
the b n and c n and clear denominators, then we obtain the relation 
/ N \ / N \ N 

E b " Tn E c « Tn_1 = E ^ rn_1 ( mod TN+l )- 

\n=l J \n=0 / n=l 

The lemma follows by expanding the left side and comparing the coefficients on each side. □ 
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As stated above, L(T,E/K) satisfies a functional equation: there is e(E/K) G {±1} such that 

L(T,E/K) = e(E/K) ■ (qT) N ■ L(l/(q 2 T),E/K), (2.5) 
hence we have the relation 

c n = e(E/K)-q 2n - N ■ c N . n , 0<n<N. (2.6) 

In the following lemma we write (^) for the Legendre symbol in of e w , defined in Table CEJ and 
give a formula for e(E/K) (cf. corollary 5 of [Hj). 

Lemma 2.3. 

e{E/K) = (-l)* M+ .\\(^-). 

Proof. This follows from calculations in [Rj where the sign is the global root number of E/K and 
is given by a product of local root numbers. If tt G \U\, then the local root number is trivial by 
proposition 8 of loc. cit. with r = 1. If tt G M U A, then we apply parts (ii) and (iii) theorem 2 of 
loc. cit. with t = 1 for the remaining cases. Note, if tt G A and does not have Kodaira symbol I* 
with n > 0, then we need the assumption that q is not divisible by 2 or 3. □ 

We observe that the recursive relation given by lemma 12.21 enables us to perform a consistency 
check when trying to compute L(T,E/K): the b m and c m are integers, so for each n > 1, the 
integer Ysm=i ^ m ' c n-m must be divisible by n. A second consistency check is to compute c n for 
one or more n > |_7V/2j using the same method as for smaller n and then to verify that (I2.6P holds. 
While one would not want to use the latter check when computing large data sets, it is very useful 
for making sure that the calculations are correct because one can use it to test a small subset of 
data. 



2.3 Quadratic Twists 

We continue the notation of the previous sections. Thus we fix an elliptic curve E/K and a minimal 
Weierstrass model y 2 = x 3 + ax + b of E over Ok- For each / G K x , we define Ef/K to be the 
elliptic curve with affine model y 2 = x 3 + f 2 ax + / 3 6. 

Lemma 2.4. Suppose L/K is an extension. If an elliptic curve over K is L-isomorphic to E, 
then it is K -isomorphic to some Ef and y/J G L. Conversely, if yfj G L, then E and Ef are 
L-isomorphic. 

Proof. If y 2 = x s + a'x + b' is an affine model for an elliptic curve E' /K, then an L-isomorphism 
E —7- E' must take the form (x,y) h-> (x/c 2 ,y/c 3 ) for some c G L x (see page 50 of [S]); recall j is 
non-constant, hence neither nor 1728. In particular, a' = c 4 a and b' = c e b, so / = c 2 = b'/a' lies 
in K x and E' = Ef. Conversely, if c = ±\/7 £ L, then (x,y) i— >■ (x/c 2 ,y/c 3 ) is an L-isomorphism 
E^t-Ef. □ 
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The lemma implies L = K(y/J) is the smallest extension over which E,Ef are L-isomorphic. Hence 
if / lies in the complement K x — (K x ) 2 , then Ef/K is a so-called quadratic twist of E/K. 

Lemma 2.5. Ef, E g are K -isomorphic if and only if f = gc 2 for some c G K x . 

Proof. Replace E/K by E g /K and apply the previous lemma with L = K. □ 

We define the family of polynomials 

3~ = {/ G Or '■ f is monic, square-free, prime to A} 

and write JjCJ for the subset of / satisfying deg(/) = d. The previous lemma implies the Ef 
are pairwise non-X-isomorphic for / 6 3", while for a fixed d we will see that deg(L(T, Ef/K)) is 
independent of / G 3^. The latter fact would not be true if we dropped the condition that / be 
relatively prime to A. As we will see, if a G ¥ x is a non-square, then L(T, E a f/K) = L(—T, Ef/K), 
hence the reason we restrict to monic /. 

If / G 3~, then one can easily verify that y 2 = x 3 + f 2 ax + f 3 b is a minimal Weierstrass model 
for Ef over 0^ with discriminant / 6 • A. There is a unique integer e such that the substitution 
(t,x,y) i->- (l/u,x/u e ,y/u e ) yields a minimal Weierstrass model for Ef over F g [ii], and we glue the 
models together over P 1 — {0, oo} to construct the minimal Weierstrass model £/ — >■ P 1 . We write 
Mf and Af respectively for the divisors of multiplicative and additive reduction respectively of 

We write A 1 for the complement P 1 — {oo}. If tt G (A 1 ), then one can easily verify that 

M f n A 1 = M n A 1 , Af n A 1 = (A n A 1 ) U {tt G |A x | : tt|/}. 

If vr G Mj n A 1 , then one can also easily verify that £ f /¥ n has the same splitting behavior as £ /¥ n 
if and only if the image of / is a square in F„-, and otherwise it has the opposite splitting behavior; 
that is, 

Mf n A 1 = |tt G M ± n A 1 : (J-} = ill U |vr G M T n A 1 : Q = . 

If / G 3~d and d is even, then £ and £/ are isomorphic over Fqo. On the other hand, if d is 
odd, then the Kodaira symbols of £/Foo and £//Foo form an unordered pair {S, S*}, where S G 
{I n , II, III, IV}. 

If we fix a non-square a G F* and / G 3", then a similar calculation for E a f shows that the Kodaira 
symbols for £j and £ a j are the same for all tt G (P 1 !, so M af = M f and A af = A f . If vr G A f , 
then £,f and £ a j are isomorphic over F^. On the other hand, for every tt G (P 1 — ^4/|, there is a 
unique quadratic twist of Ef/¥ w , which we call the scalar twist, and it is easy to show that E a f/¥ v 
is isomorphic to £//F T if deg(7r) = [F,r : ¥ q ] is even, otherwise it is the scalar twist. We call E a f 
the scalar twist of Ef. 
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2.4 L-functions of Quadratic Twists 

We continue the notation of the previous section and fix an elliptic curve E/K and a quadratic twist 
Ef/K. If Uf, Mf, and Af are the primes over which Ef has good, multiplicative, and additive 
reduction respectively, then we can use the results of section 12.21 to infer that the L-function 
L(T,Ef/K) has Euler product expansion 

L(T,E f /K)= 11 L^^^f/W^Y 1 ■ Y[ fl-T^W)" 1 - ]J (l + T deg W)~\ 

There is an important observation which relates this to the Euler product expansion in (|2.2[) of 
L(T,E/K): if tt lies in U D Uf and if Xn(f) G {if} denotes the Legendre symbol of / (mod 7r), 
then 

L(r,£ / /F 7r ) = L(x 7r (/)T,£/F 7r ). 

In particular, if one has precomputed sufficiently many Euler factors for L(T, E/K), then for most 
of the Euler factors of L(T,Ef/K), the cost of computing the factor is essentially the cost of 
computing xM)- 

2.5 Pullbacks 

We continue the notation of previous sections and fix an elliptic curve E/K and a minimal Weier- 
strass model y 2 = x 3 + ax + b of E over Ok- We write M, A respectively for the subsets of primes 
in K over which E has multiplicative and additive reduction respectively. 

We fix a rational function field L = ¥ q (w) and a non-constant element 9 G L, and we write 
9* : K —7- L for the embedding induced by sending t to 9. Let Ol C L be the integral closure of 
Ok, and let Oob C L be the integral closure of Ooo C K. We call the primes of Ol the finite primes 
of L and the primes of Ooo the infinite primes of L. In general, Ol is not Fgfu;] and Oqo is not the 
local ring with uniformizer 1/w, but rather the infinite primes are the poles of 9. If tt is a prime 
(finite or infinite) in L, we write 9(ir) for the corresponding prime in K, for the inertial degree 
of tt over 9(tt), and for the ramification degree. 

We write El for the elliptic curve over L (eliding the dependence on the choice of 9) with affine 
model y 2 = x 3 + 9*{a)x + 9*(b), and thus the coefficients of the model are rational functions in w 
when a, b are viewed as elements of L via 9*. A priori, this model is not a minimal Weierstrass 
model of El over Ol, but if tt is a finite prime that is unramified in L or if 9{tt) does not lie in 
A, then the model is a minimal Weierstrass model over the local ring O^. Similarly, if tt £ Oob is 
an infinite prime, then a minimal Weierstrass model for E over 0^ is guaranteed to be a minimal 
Weierstrass model over 0^ only if tt is unramified over 9(tt) or if oo does not lie in A. 

Now suppose that tt is a prime of L such that tt is ramified over 9(tt) and -E has bad reduction 
over #(7r), and let y 2 = x 3 + a^x + b^ be a minimal Weierstrass model of E 1 over 0^) and let 
Ag^) be the discriminant of this model. If E has Kodaira type I n over 9{tt) and if e = (e 
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being the unique integer used to obtain a model for E over F g [u]), then y 2 = x 3 + 6*(a 1T )x + 9*{b v ) 
is a minimal Weierstrass model of El over 0^ with discriminant A n = 9*(Ag^) and El has 
Kodaira type I en over n. On the other hand, if 0(ir) lies in A, then the Kodaira type of El over 
7T may differ from the Kodaira type of E over 6(ir) depending on e = e n . More precisely, if E 
has Kodaira type I* over 9(tt), then then El has Kodaira type I en or I* n over 7r if e is even or 
odd respectively. Otherwise, the discriminant Ajr for a minimal Weierstrass model of El over 0^ 
satisfies ord 7r (A 7r ) = e-orde( 7r )(Ae( 7r )) (mod 12), hence the Kodaira type of El over ir is completely 
determined by the Kodaira type of E over 8(ir) and table [TJ 

Aside from the fact that one can use pullbacks to generate new elliptic surfaces from old, the other 
important role they play lies in the fact that any elliptic curve over L with non-constant j-invariant 
can be written as a quadratic twist of the pullback of the 'versal' elliptic curve E/K with affine 
model 

2 o 108t 432t 
y = x x + 



t - 1728 t - 1728 

One can easily verify that this elliptic curve has j-invariant t, hence for an elliptic curve over L one 
can take 9 to be the j-invariant and use lemma [231 to infer that an appropriate quadratic twist of 
the pullback will be the original elliptic curve over L. We remark that if y/2 F q , this model is the 
twist of the model in [Si proof of prop. III. 1.4] by the quadratic extension K(s/2)/K. The latter 
model has split-multiplicative reduction at t = oo while our model forces e(E/K) = —1. In both 
cases the L-function has degree one and thus in our model we have L(T, E/K) = 1 — qT. One can 
verify that the point P = (4,8) lies on E and has height 1/2, thus the Mordell-Weil and analytic 
ranks of E are both one. 



2.6 L-functions of Pullbacks 

We continue the notation of the previous section and fix an elliptic curve E/K and a pullback 
El/L. If Ul, Ml, and Al are the primes over which El has good, multiplicative, and additive 
reduction respectively, then the L-function L(T, El/L) has Euler product expansion 

L(T,E L /L)= 11 LiT^W^L/Wv)- 1 ■ Yl (1 - T**^)- 1 • (l + T^W) -1 . 

As in the case of quadratic twists, if one has computed enough information for E/K, then it is 
relatively cheap to compute most of the Euler factors of L(T,El/L). More precisely, if El has 
good reduction over ir and if E has good reduction over p = 6(tt), then 

L(T, £ L /F,) = 1 - V T + g dc sMr 2 , 

where f n is the inertia degree of ir over p and a p / w can be determined by the expansion (12.ip . In 
practice, it is easier to keep track of a p ™ for several n rather than use (|2.ip directly because, among 
other reasons, it is difficult to compare elements of F p with the corresponding subfield of F,,- = F p / W . 
Nonetheless, the additional work one must do is small since 'most' elements of F q n have degree n 
over F q , and the upshot is that most of the work of computing b n in the corresponding expansion 
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(1231) for L(T,E L /L) is the cost of explicitly evaluating the map 6 : F 1 (F 9 n) -+ P 1 ^) for all 
elements in the domain. 



3 ELLFF 

The discussion in section 12.21 above naturally leads to a naive algorithm to compute the L- function 
of a non-isotrivial elliptic curve defined over a function field via counting points on a finite number 
of its fibers. Moreover, if sufficiently many Euler factors have been computed for the versal elliptic 
curve, one can realize a given elliptic curve as a pullback and quadratic twist and use the results of 
section [2T31 - l2.6l to significantly reduce the number of fibers on which one needs to count points. The 
authors have written a standalone C++ library called ELLFF built upon Shoup's Number Theory 
Library [NTLJ, which can be added as a module to the free open-source mathematics software 



system Sage S + 10 . The package allows anyone to efficiently compute these L-functions on their 
own|§ 

Internally, the library uses tables, computed on demand, to represent Euler factors. If one asks for 
the L-function of a curve, then the library demands the minimal number of tables necessary. One 
may also demand and manipulate the tables directly, e.g. in order to study how the sizes of special 
fibers vary. The library uses a database in order to reduce the complexity of computing a table, 
e.g. by returning a previously calculated copy of the table or twisting a table for another curve with 
the same j-invariant. A user whose database has the appropriate tables for the versal curve will 
benefit from such reductions, and thus we have made available a modest collection for download. 
The user may also save their own tables in the database in order to facilitate calculating tables for 
families of curves. For more information and setup instructions, see http : / / ellf f . sagemath . org. 



4 Computations 

Using the discussion from section [2j a database of L-functions was amassed for the family of 
quadratic twists of the following four elliptic curves (with notation consistent with that found in 
pP]): 

^222 :y 2 = x 3 - 27(t 4 - t 3 + t 2 )x + 27(2t 6 - 3t 5 - 3t 4 + 2t 3 ), 

X 211 :y 2 = x 3 - 27t 4 x + 54t 5 (t - 2), 

X 321 :y 2 = x 3 - W8t 3 (4t - 3)x + 432t 5 (8t - 9), 

X 43 i : y 2 = x 3 - t 3 (27t - 2A)x + t 4 (54t 2 - 72t + 16). 

These are the only elliptic curves, up to isogeny, over ¥ q (t) with (q,6) = 1 such that #M = 2 and 
#A = 1. They are normalized so that 00 G M + , t £ A, and t — 1 S M, forcing the L-function 

2 The library currently only allows for characteristic not 2 or 3, though handling these cases is straightforward and 
will be addressed in a future release. 
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to be trivial for each of the curves. Note that the first curve is the Legendre curvejfl given by the 
alternative model 

A 22 2 : y 2 = x(x + t)(x + t 2 ). 

For each of these curves, we considered all prime q G Q = {5,7,... ,29} and computed the L- 
functions of all the twists with bounded degree. The bound on the degree was determined by 
considerations on computational feasibility and depended on the size of the field of constants. The 
following table lists the number of twists over ¥ q of degree d in each family 3"^. This number does 
not depend on which of the four curves above one considers. A blank entry in the table denotes 
that the L-functions for all twists for the given d and q were not determined due to the computation 
requiring an excessive amount of time. 





5 


7 


11 


13 


17 


19 


23 


29 




3 


5 


9 


11 


15 


17 


21 


27 




13 


31 


91 


133 


241 


307 


463 


757 




71 


227 


1019 


1751 


4127 


5867 


10691 


22007 




345 


1573 


11181 


22729 


70113 


111421 


213762 


638121 




1739 


11033 


123029 


295523 












8677 


77203 
















43407 


















217009 


















1085075 
















All 


1356339 


90072 


135329 


320147 


74496 


117612 


224937 


660912 



Table 2: Number of Twists in 3^ for q S Q 



4.1 Goldfeld's Conjecture 

In 1979 Goldfeld [G] conjectured an average value for the analytic rank of a family of quadratic 
twists of a fixed elliptic curve 



Conjecture 4.1 (Goldfeld). For D a discriminant, 

Km Z»<Drm = 1 

D->oc #{d : \d\ < D} 2 y J 

where r(Ed) is the order of vanishing at s = 1 of the L-function of the quadratic twist E^/Q- 



Strictly speaking, X222 is a twist by —t of the usual Legendre curve model y = x(x — l)(x — t). 
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Goldfeld's conjecture concerns the analytic rank of an elliptic curve, though it is important to 
note that many authors replace the analytic rank with the algebraic rank (i.e. the rank, as a free 
Z-module, of the group E(K) of i^-rational points on E modulo torsion), invoking the Birch and 
Swinnerton-Dyer conjecture if needed. For a survey of results on the average value and variation 
of the (algebraic) ranks of elliptic curves in a family of quadratic twists in the number field setting, 
sec [RSI. A more recent paper [BMSW] provides data for the average value and distribution of 
the analytic ranks of elliptic curves over Q ordered by conductor. Thus the reader should be wary 
of concluding that the data presented therein either supports or undermines Goldfeld's conjecture, 
which considers the family of quadratic twists of a fixed elliptic curve and not all elliptic curves 
with bounded conductor. 

Goldfeld's conjecture has a direct analog in the function field setting: for an elliptic curve E over 
K, we set its analytic rank r to be the order of vanishing of L(T, E/K) at T = 1/q. Instead of 
considering all twists by d with \d\ < D, we consider those twists in 3* D = {J d<D 3~d and let D grow 
to infinity as before: 

Conjecture 4.2. For D a positive number, 

where r(Ef) is the order of vanishing at s = 1 of the L-function of the quadratic twist Ef/¥ q (t). 

One would like a lower bound on the average analytic rank over the family of interest 3^ analogous 
to that found in Proposition 1 on p. 114 of [G]. Unlike this proposition where the average is taken 
over all discriminants, determining the average over 3^ is non-trivial. But using the functional 
equation, it is clear that if in the limit the average of the sign of the functional equation over 3^ is 
0, the average analytic rank over 3^ is at least 1/2. We next prove such a lower bound using this 
line of argument. 

We begin by letting M n A 1 = {711, . . . 7r r } be the finite primes where E/K has multiplicative 
reduction and setting N = 7Ti • • -ir r . 

Proposition 4.3. There exists Ed 6 {±1} such that for all f G 3"^, 

e(E f /K) = e d -e(E/K)-(Jj\, (4.3) 
where (-^f) is the Jacobi symbol of N. 

Proof. We proceed by examining the contribution to the sign from the places of bad reduction. 

Case 1: If 7r £ A\ n A 1 , E/K and Ef /K have the same Kodaira type at tt. Thus there is no change 
to the local contribution from e(E/K) to e(Ef/K) for such ir. 
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Case 2: If tt is a finite prime that divides /, then Ef has type Iq reduction over tt. Thus the 
contribution to the sign in this case is given by 

'^pj = (ZL\ = q ^ ( m od 4), 

implying the total contribution e/ to the sign coming from those tt G Af — A satisfies 

ef = q d (mod 4). 

Thus the change in the local contribution from e(E/K) to e(Ef/K) from these primes depends 
only on d. 

Case 3: For tt E Mf n A 1 = M n A 1 , the splitness at tt changes if and only if ^ j = —1. Thus the 
total change in the local contribution from e(E/K) to e(Ef/K) is (4f 

Case 4-' For tt = oo, the reduction of E/K and Ef/K are the same if d is even by the discussion 
in section 12.31 If d is odd, then the reduction depends only the leading coefficient of /. Thus the 
local contribution to the sign is independent of / G 9^, so the change in the local contribution from 
e(E/K) to e(Ef/K) for tt = oo depends only on d. 

These cases exhaust all possible changes to the sign of the functional equation introduced by 
twisting, yielding equation 14,31 □ 

Corollary 4.4. 

Corollary 14.41 reduces the average of the sign of the functional equation to the average of the Jacobi 
symbol over 3^- The following proposition is due to private correspondence with E. Kowalski: 



Proposition 4.5 (Kowalski). With notation as above, we have 

E/e7 d (jf 
hm — - 

d-+oc #Jd 

Proof. Unless stated otherwise, we write f,g,h G ¥ q [t] for arbitrary monic polynomials. Write 
A = NN' and let xa(/) be the characteristic function of those / that are square- free and coprime 
to A. Setting 



we then have 



deg(/)=d V 7 deg(/)=d 
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Let jtx(-) be the Mobius function for polynomials. If deg(g) > 0, then = 0> an( i otherwise 

^Zh\ g ^{h) = 1. Thus / i-> X^3 2 |/ A*(f) an d / ^ Sh|(A /) mC 1 ) are the characteristic functions for 
square-free polynomials and polynomials coprime with A respectively, and hence 

XA(/) = X>fo) E MW- 

9 2 l/ fc|(A,/) 

Note, if (g, A) 7^ 1, then the right sum over h vanishes, hence we can restrict to g such that 
(g, A) = 1. In particular, if we substitute into the above expression for Ad and rearrange terms we 
have 

A d = em/o e e (T)- 

h\A dcg( 9 )<§, (<?,A)=1 9 2 /*l/, deg(/)=d V 7 

If we write / = f\g 2 h in the innermost sum, then we have 

7i 



4* = E E ^) E 

h|A V 7 deg( ff )<f, (s,A)=l dcg(/i)=d-2dcg( 9 )-dcg(fe) 



A/ 



Moreover, if we write B e for the sum B e = ^deg(/)=e (iv) anc ^ ^ we suppose e > deg(iV), then 

«- - E (£) E 1 = E (£) " e - d ' 8, " ) = ° 

a€F,[t]/(JV) dcg(/) = e a£F,[t]/(JV) 
/ = a(mod N) 

(because the last sum is a complete character sum). Therefore if we write e = d — 25 — deg(/i) and 
suppose e < deg(iV) we have 

A2 = E^)( a t) £ £ Ka) B d _ 2S _ dcg(h) . 

h\A V 7 i( d _ d eg(jV)-deg(fc))<<5<§ dcg( 9 )=<5, ( 9 ,A)=1 

Observe that for any e, 5 > 0, we have |5 e | < g e and (^2deg(g)=5 1) — 9 5 > thus for d > 1 

lAfi < E^) £ g*- 2 *-***) e 1 

ft|A i( ( i_dcg(AT)-dcg(h))<<5<f deg( 9 )=5, (<?,A)=1 

< e E ^ deg(JV) ? d/2 

fe|A ±(<i-deg(iV)-deg(/i))<<5<§ 



< q dcg(N)+d/2 

where the implied constant depends on A and N. On the other hand, #3"^ » q d , so 



proving the propsition. □ 



Jim iM « lim ^(iv)-^ = 0> 



This proposition then leads to the desired corollary: 
Corollary 4.6. With notation as above, we have 
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4.2 Average Analytic Rank Data 



We define 



E 



iE f ) 



D 



to be the average rank of the family of quadratic twists of E up to degree D. This value was 
calculated for the four elliptic curves discussed above with increasing D, and the data is presented 
in tables [SHE] below, where the dependence of the average rank on D is made explicit. As in the 
case of Table an empty entry denotes that those computations were not done. 



D 



7 



11 



13 



17 



19 23 



2!) 



1 1.000 0.400 0.667 0.636 0.733 0.588 0.571 0.704 

2 0.688 0.667 0.680 0.674 0.668 0.679 0.661 0.652 

3 0.644 0.669 0.622 0.629 0.610 0.607 0.588 0.576 

4 0.653 0.659 0.638 0.620 0.605 0.599 0.588 0.575 

5 0.666 0.633 0.590 0.581 

6 0.628 0.609 

7 0.623 

8 0.592 

9 0.582 

Table 3: n(X 222 ,D) for q G Q 



D 



7 



11 



13 17 19 23 



2!) 



1 0.333 0.400 0.444 0.636 0.467 0.588 0.571 0.481 

2 0.688 0.556 0.540 0.549 0.527 0.568 0.562 0.545 

3 0.598 0.601 0.533 0.579 0.536 0.562 0.526 0.524 

4 0.662 0.562 0.543 0.534 0.529 0.532 0.525 0.521 

5 0.586 0.565 0.525 0.538 

6 0.634 0.539 

7 0.554 

8 0.581 

9 0.535 

Table 4: fi(X 2n ,D) for q G Q 



D 



7 



11 



13 17 19 23 



2!) 



D 



11 



13 



17 19 23 29 



1 0.333 0.400 0.444 0.636 0.600 0.471 0.571 0.481 

2 0.562 0.556 0.540 0.618 0.590 0.580 0.587 0.585 

3 0.690 0.570 0.577 0.605 0.583 0.565 0.552 0.558 

4 0.625 0.609 0.569 0.574 0.559 0.555 0.551 0.543 

5 0.618 0.571 0.553 0.554 

6 0.602 0.569 

7 0.587 

8 0.568 

9 0.556 

Table 5: ^(X 321 ,D) for q G Q 



1 0.333 0.800 0.444 0.636 0.467 0.588 0.571 0.556 

2 0.562 0.611 0.640 0.618 0.613 0.611 0.616 0.614 

3 0.621 0.646 0.602 0.632 0.580 0.594 0.567 0.557 

4 0.616 0.617 0.593 0.590 0.576 0.575 0.562 0.553 

5 0.592 0.600 0.558 0.555 

6 0.601 0.574 

7 0.575 

8 0.568 

9 0.548 

Table 6: fi(X 431 ,D) for q G Q 



Considering each table separately, the individual columns present the data pertaining to Goldfeld's 
conjecture. In particular, for the largest data sets with q = 5, there is a slow convergence to 
the conjecture value of 1/2. On the other hand, each row of a table presents data relevant to 
Katz-Sarnak [KS], where one lets q grow to infinity to determine that the conjugacy classes of the 
Frobenius automorphism are equidistributed in the special orthogonal group of N x N matrices 
with respect to Haar measure, where N is the degree of the L-function. 
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We can also consider how the average ranks varies between each of the four curves for a fixed q, 
as presented in figures [1] through [8] below. Recall that the four cures are not isogenous but have 
nearly the same reduction types. Even for the smallest data sets (q £ {17, 19, 23, 29}), there is good 
numerical evidence that the average ranks for each of the four curves are converging to the same 
value for any given q. Again, q = 5 provides the strongest evidence that this value is 1/2. Note 
also that in general the average rank of the Legendre curve X222 dominates the rank of the other 
three curves. 




0.3 - 



0.2 - 



0.1 - 



I i i i i i i i I 

1 23456789 

D 

Figure 1: Variation of fi(X{, D) as i varies for q = 5 
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17 




18 
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4.3 Rank Distributions 



Combining Goldfeld's conjecture with the parity conjecture leads to a conjecture on the density of 
ranks in a family of quadratic twists (for details of the formulation, see section 7.6 in [RSJ): 

Conjecture 4.7. With notation as in Conjecture \4-S\ 

W G *b ■ r ( E f) = °} _ ^ #{/ e ff£ : r(tf,) = 1} _ 1 



whereas 



lim — = lim 

D^oo #3% D^oo #3* 2 



lim #Ve*b:r(E,)>2} =0 _ 



As the data on average rank above suggests, the distribution of analytic ranks for our four families 
is close to that predicted by the density conjecture with a non-trivial number of twists with rank 
greater than or equal to twoQ We present the relevant data in the table below, where we have 
removed the dependence of the distribution on the degree of the twisting polynomials and instead 
consider all the L-functions we were able to compute given some q. For the dependence of the rank 
distribution on the degree, see the tables in the appendix. 



4 The largest rank discovered was a rank 5 curve, a twist of X222 /F 5 (£) by / = t 7 + 2t 6 + t 5 + 4t 4 + At 3 + 1 2 + 2t + 1. 
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Q 

Rank 


5 

1 2 > 3 


7 

1 2 > 3 


11 

1 2 > 3 


13 

1 2 > 3 


A222 
-^211 
-^321 


.461 .498 .039 .002 
.483 .500 .018 .000 
.473 .499 .027 .001 
.477 .499 .023 .001 


.447 .499 .053 .002 
.481 .500 .019 .001 
.468 .497 .031 .003 
.464 .499 .036 .001 


.457 .498 .043 .002 
.488 .500 .012 .000 
.474 .500 .026 .000 
.471 .500 .029 .000 


.462 .498 .038 .002 
.481 .500 .019 .000 
.474 .499 .026 .000 
.474 .499 .026 .001 




Q 

Rank 


17 

1 2 > 3 


19 

1 2 > 3 


23 

1 2 > 3 


29 

1 2 > 3 


-^222 
A^211 
A321 
X431 


.450 .498 .059 .002 
.485 .500 .014 .000 
.471 .500 .029 .000 
.463 .499 .037 .000 


.452 .498 .048 .002 
.484 .500 .016 .000 
.473 .500 .027 .001 
.463 .499 .036 .001 


.458 .498 .042 .002 
.488 .500 .012 .000 
.475 .500 .025 .000 
.470 .499 .030 .001 


.463 .499 .036 .001 
.490 .500 .010 .000 
.479 .500 .021 .000 
.474 .499 .026 .001 



Table 7: Rank Distributions for All Curves over All d and q £ Q 



5 Conclusion 

The remarkable property that the L-function of a non-isotrivial elliptic curve over a function field 
is a polynomial yields an effective algorithm to determine its coefficients by computing the number 
of points on a finite number of fibers. These fibers precisely correspond to the Euler factors that 
determine L(E/K,T), and by realizing a given curve as a quadratic twist or pullback of another 
curve, the number of Euler factors that need to be computed can be minimized. In particular, 
the versal elliptic curve provides a (non-canonical) choice for an elliptic curve from which one can 
pullback and twist to recover any given elliptic curve, allowing for the efficient computation of the 
given curve's L-function provided sufficiently many Euler factors have been precomputed. These 
algorithms have been developed into ELLFF, a software library for the open-source mathematical 
software system Sage, allowing anyone to quickly compute such L-functions. 

Experimentally, we computed the L-functions of four different families of quadratic twists in order 
to examine their analytic ranks for numerical evidence pertaining to Goldfeld's conjecture. Using an 
elementary argument, we know the asymptotic average rank over our family of quadratic twists is 
at least 1/2 as the degree of the twisting polynomial becomes arbitrarily large. Unlike the situation 
in number fields, the case of function fields provides strong evidence, especially for the largest 
data sets, that this average is indeed 1/2, thus supporting the validity of Goldfeld's conjecture in 
the function field case. Moreover, the experimental data also suggests that the analytic ranks are 
distributed closely with the density conjecture's prediction. Nonetheless the presence of a non- 
trivial amount of curves of rank at least 2 in even the largest data sets may suggest the convergence 
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to this distribution is rather slow. 



This work is part of a small but growing body of computational number theory directly focused on 
function fields. Historically, computational number theorists have primarily worked over number 
fields, in particular Q. This (understandable) bias has produced a dearth of algorithms and data 
for the function field setting, despite the fact that many of the ideas from number fields can be 
formulated more generally for any global field. There is much work left to be done - both theoretical 
and computational - for the case of function fields, but we believe the example of L-functions of 
elliptic curves indicates that the effort is worthwhile and yields interesting mathematics. 
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Appendix: Analytic Rank Distribution Tables 

The following tables give the distribution of analytic ranks, making explicit their dependence on 
the degree d of twists considered in 9^. 







Rank 


o 


1 


2 


3 


4 


5 


Q = 


5 


?i 


0.333333 


0.333333 


0.333333 














0.384615 


0.615385 
















0.422535 


0.521127 


0.056338 












4 


0.417391 


n 510145 


079464 












J 5 


0.420357 


f! 4Q6969 


077631 












To 
J 6 


0.444163 


4Q79Q9 


n 054Q73 


003349 


000930 








J 7 


449509 


n 49701 7 


n 0571 57 


003087 


000907 


000093 






3*8 


0.459645 


0.498025 


0.040256 


0.002009 


0.000065 








3"9 


0.461748 


0.498378 


0.038150 


0.001628 


0.000095 








All 


0.460615 


0.498274 


0.039264 


0.001751 


0.000094 


0.000001 


<? — 


; 7 


J 1 


0.600000 


0.400000 














To 


0.354839 


0.580645 


0.064516 












To 
J 3 


0.418502 


0.502203 


0.070485 


0.008811 










- 1 -1 


0.419580 


0.503497 


0.076923 












T 5 


0.438684 


0.497326 


0.060818 


0.003172 










T 6 


0.448610 


0.498685 


0.051280 


0.001425 










All 


0.446787 


0.498634 


0.052947 


0.001632 






Q — 




J 1 


0.444444 


444444 


0.111111 












J 2 


0.384615 


0.549451 


0.065934 












T 3 


0.442591 


0.501472 


0.052993 


0.002944 










T 4 


0.433056 


0.497809 


0.066005 


0.003041 


0.000089 








T 5 


0.459054 


0.498232 


0.040828 


0.001845 


0.000041 








All 


0.456731 


0.498252 


0.043021 


0.001951 


0.000044 




Q — 


13 


J 1 


0.454545 


0.454545 


0.090909 












T 2 


0.398496 


0.533835 


0.060150 


0.007519 










3*3 


0.439749 


0.499143 


0.057110 


0.003998 










3*4 


0.442430 


0.497998 


0.057064 


0.002508 












0.463263 


0.498073 


0.036630 


0.001966 


0.000068 








All 


0.461629 


0.498087 


0.038204 


0.002018 


0.000062 






17 


3*1 


0.400000 


0.466667 


0.133333 












3*2 


0.406639 


0.526971 


0.062241 


0.004149 










3" 3 


0.452144 


0.495517 


0.046038 


0.006300 










T 4 


0.449574 


0.498210 


0.050162 


0.002011 


0.000043 








All 


0.449568 


0.498148 


0.049989 


0.002255 


0.000040 




<? — 


19 


3*1 


0.470588 


0.470588 


0.058824 












3-2 


0.400651 


0.521173 


0.071661 


0.006515 










3 r 3 


0.448440 


0.500767 


0.050111 


0.000682 










T 4 


0.452392 


0.498236 


0.047451 


0.001921 










All 


0.452063 


0.498419 


0.047648 


0.001871 






<? = 


23 




0.476190 


0.476190 


0.047619 












3-2 


0.406048 


0.522678 


0.071274 












*3 


0.458984 


0.499111 


0.039940 


0.001871 


0.000094 








T 4 


0.457616 


0.498264 


0.042271 


0.001824 


0.000023 








All 


0.457577 


0.498353 


0.042221 


0.001823 


0.000027 




Q = 


29 


3-1 


0.407407 


0.481481 


0.111111 












3"2 


0.418758 


0.515192 


0.063408 


0.002642 










3 3 


0.465170 


0.498705 


0.034216 


0.001908 










T 4 


0.463392 


0.498916 


0.036535 


0.001127 


0.000030 








All 


0.463398 


0.498927 


0.036492 


0.001154 


0.000029 





Table 8: Rank Distribution of X 22 2/%(t) 
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Rank 12 3 4 

9 = 5 31 0.666667 0.333333 

31 0.307692 0.615385 0.076923 

31 0.450704 0.521127 0.028169 

31 0.405797 0.510145 0.084058 

31 0.465210 0.502013 0.032777 

3" 6 0.427337 0.499942 0.072030 0.000691 

3> 0.483125 0.499758 0.016748 0.000369 

3" 8 0.456834 0.499611 0.043132 0.000424 

3" 9 0.488469 0.499751 0.011520 0.000256 0.000005 

All 0.482791 0.499737 0.017179 0.000289 0.000004 

9 = 7 31 0.600000 0.400000 

31 0.419355 0.580645 

3-3 0.440529 0.511013 0.048458 

31 0.471710 0.502225 0.024793 0.001271 

3% 0.467597 0.499864 0.031904 0.000634 

3" 6 0.482935 0.499566 0.016955 0.000544 

All 0.480738 0.499700 0.018996 0.000566 

q = 11 31 0.555556 0.444444 

3- 2 0.450549 0.549451 

31 0.481845 0.504416 0.013739 

31 0.479116 0.499419 0.020034 0.001431 

3% 0.488324 0.499963 0.011599 0.000114 

All 0.487493 0.499982 0.012303 0.000222 

q = 13 31 0.454545 0.454545 0.090909 

3- 2 0.458647 0.541353 

31 0.457453 0.503141 0.039406 

31 0.484667 0.500506 0.014827 

3% 0.480937 0.499809 0.019024 0.000230 

All 0.481063 0.499892 0.018832 0.000212 

q = 17 31 0.533333 0.466667 

31 0.468880 0.531120 

31 0.480979 0.501817 0.017204 

31 0.485716 0.500221 0.014063 

All 0.485409 0.500403 0.014189 

q = 19 31 0.470588 0.470588 0.058824 

31 0.452769 0.527687 0.019544 

3-3 0.468723 0.501108 0.029828 0.000341 

3-4 0.484918 0.500157 0.014925 

All 0.484024 0.500272 0.015687 0.000017 

q = 23 31 0.476190 0.476190 0.047619 

31 0.457883 0.522678 0.019438 

31 0.487513 0.500889 0.011505 0.000094 

31 0.487659 0.500089 0.012252 

All 0.487590 0.500171 0.012235 0.000004 

q = 29 31 0.518519 0.481481 

31 0.467635 0.517834 0.014531 

3-3 0.488345 0.500523 0.011042 0.000091 

31 0.489674 0.499968 0.010283 0.000075 

All 0.489605 0.500006 0.010313 0.000076 



Table 9: Rank Distribution of X 2 n/¥ q {t) 



24 



Rank 12 3 4 

9 = 5 31 0.666667 0.333333 

?2 0.384615 0.615385 

3 3 0.380282 0.521127 0.098592 

31 0.446377 0.504348 0.043478 0.005797 

3% 0.444508 0.498562 0.053479 0.003450 

y 6 0.452345 0.498790 0.047021 0.001844 

0*7 0.460525 0.497892 0.039302 0.002235 0.000046 

3* 8 0.469898 0.498468 0.030059 0.001567 0.000009 

3*g 0.474192 0.499060 0.025765 0.000947 0.000036 

All 0.472877 0.498928 0.027065 0.001098 0.000032 

9 = 7 31 0.600000 0.400000 

3*2 0.419355 0.580645 

3*3 0.458150 0.511013 0.030837 

3* 4 0.448188 0.495868 0.048315 0.007629 

3% 0.468504 0.499592 0.030998 0.000906 

3*6 0.0.468881 0.496872 0.030970 0.003238 0.000039 

All 0.468436 0.497247 0.031264 0.003020 0.000033 

q = 11 3"! 0.555556 0.444444 

3- 2 0.450549 0.549451 

3*3 0.457311 0.504416 0.038273 

3*4 0.465254 0.500850 0.033897 

3% 0.475002 0.499614 0.024921 0.000463 

All 0.474052 0.499782 0.025745 0.000421 

q = 13 31 0.454545 0.454545 0.090909 

s 3" 2 0.421053 0.541353 0.037594 

3*3 0.448886 0.500857 0.047973 0.002284 

3* 4 0.464297 0.500506 0.035197 

3*5 0.0.474582 0.499125 0.025379 0.000914 

All 0.473689 0.499249 0.026207 0.000856 

q = 17 31 0.466667 0.466667 0.066667 

3- 2 0.439834 0.531120 0.029046 

3*3 0.459414 0.500363 0.038769 0.001454 

31 0.471696 0.499594 0.028083 0.000628 

All 0.470911 0.499732 0.028686 0.000671 

q = 19 31 0.529412 0.470588 

3- 2 0.442997 0.527687 0.029316 

3- 3 0.467701 0.501108 0.030851 0.000341 

3-4 0.473385 0.499574 0.026458 0.000583 

All 0.473030 0.499719 0.026681 0.000570 

q = 23 31 0.476190 0.476190 0.047619 

3*2 0.444924 0.522678 0.032397 

3*3 0.474418 0.500889 0.024600 0.000094 

31 0.474963 0.499584 0.024939 0.000505 0.000009 

All 0.474875 0.499691 0.024940 0.000485 0.000009 

q = 29 31 0.518519 0.481481 

3- 2 0.446499 0.517834 0.035667 

3*3 0.471986 0.499659 0.027400 0.000954 

31 0.478870 0.499689 0.021082 0.000354 0.000005 

All 0.478605 0.499708 0.021308 0.000374 0.000005 



Table 10: Rank Distribution of X 321 /F q (t) 
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Rank 12 3 4 

9 = 5 31 0.666667 0.333333 

3^2 0.384615 0.615385 

3- 3 0.422535 0.521127 0.056338 

a~4 0.449275 0.498551 0.040580 0.011594 

3-5 0.456009 0.502013 0.041978 

0*8 0.452345 0.496254 0.047021 0.004379 

3*7 0.466169 0.499551 0.033704 0.000576 

3*8 0.469543 0.497307 0.030377 0.002728 0.000046 

3*9 0.478878 0.499746 0.021099 0.000261 0.000016 

All 0.476768 0.499332 0.023186 0.000695 0.000020 

9 = 7 31 0.400000 0.400000 0.200000 

3* 2 0.419355 0.580645 

3*3 0.418502 0.511013 0.070485 

3-4 0.442467 0.503497 0.054037 

3% 0.454274 0.497507 0.045228 0.002991 

3*6 0.465992 0.499268 0.033768 0.000842 0.000130 

All 0.464007 0.499178 0.035616 0.001088 0.000111 

q = 11 31 0.555556 0.444444 

3- 2 0.395604 0.549451 0.054945 

3* 3 0.448479 0.504416 0.047105 

31 0.453537 0.500850 0.045613 

3% 0.473027 0.499858 0.026896 0.000219 

All 0.471185 0.500004 0.028612 0.000200 

q = 13 31 0.454545 0.454545 0.090909 

3- 2 0.421053 0.541353 0.037594 

3*3 0.439178 0.495717 0.057681 0.007424 

31 0.457609 0.499186 0.041885 0.001320 

3*5 0.475320 0.498882 0.024611 0.001157 0.000030 

All 0.473842 0.498902 0.026026 0.001203 0.000028 

q = 17 31 0.533333 0.466667 

3- 2 0.423237 0.531120 0.045643 

3*3 0.460383 0.501333 0.037800 0.000485 

3*4 0.462981 0.499280 0.036798 0.000941 

All 0.462723 0.499490 0.036874 0.000913 

q = 19 31 0.470588 0.470588 0.058824 

3- 2 0.433225 0.524430 0.039088 0.003257 

3* 3 0.456792 0.497529 0.041759 0.003920 

3-4 0.463782 0.499260 0.036043 0.000897 0.000018 

All 0.463354 0.499235 0.036340 0.001054 0.000017 

q = 23 31 0.476190 0.476190 0.047619 

3- 2 0.429806 0.522678 0.047516 

3*3 0.467590 0.500514 0.031428 0.000468 

31 0.469808 0.499289 0.030094 0.000800 0.000009 

All 0.469620 0.499393 0.030195 0.000782 0.000009 

q = 29 31 0.481481 0.481481 0.037037 

3* 2 0.433289 0.517834 0.048877 

3-3 0.472713 0.500250 0.026673 0.000364 

31 0.474451 0.499316 0.025503 0.000727 0.000003 

All 0.474346 0.499368 0.025569 0.000714 0.000003 



Table 11: Rank Distribution of X 431 /¥ q (t) 
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